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Abstract. We report on observations of the Large Magellanic Cloud with 
the Fermi Gamma-Ray Space Telescope. The LMC is clearly detected with the 
Large Area Telescope (LAT) and for the first time the emission is spatially well 
resolved in gamma-rays. Our observations reveal that the bulk of the gamma-ray 
emission arises from the 30 Doradus region. We discuss this result in light of the 
massive star populations that are hosted in this area and address implications 
for cosmic-ray physics. We conclude by exploring the scientific potential of the 
ongoing Fermi observations on the study of high-energy phenomena in massive 
stars. 



1. Introduction 

Since the early days of high-energy gamma-ray astronomy it has become clear 
that the gamr na-ray flux receiv ed at Earth is dominated by emission from the 
Galactic disk ( Clark et al. 19681 ). This emission can be well understoo d in terms 



of cosmic-ray interactions with the interstellar medium ( Strong 20071 ). At ener- 



gies ^ 100 MeV, the generation of diffuse gamma-ray emission is dominated by 
the decay of ir^ produced in collisions between cosmic-ray nuclei and interstellar 
medium nuclei. Ultimately, the study of this hadronic gamma-ray emission may 
provide hints on the still mysterious origin of the galactic cosmic-rays. However, 
the interpretation of the galactic diffuse gamma-ray emission is complicated by 
the fact that a large number and variety of individual sources contribute along 
the line of sight to the observed emission, thus blurring the link between indi- 
vidual cosmic-ray acceleration sites and observed gamma-ray signatures in our 
Galaxy. 

Gamma rays from cosmic-ray interactions are also expected from nearby 
galaxies, and indeed, the EGRET telescope aboard the Compton Gamma-Ray 
Observatory (CGRO) has for the fi rst time detected gamm a-ray emission from 
the Large Magellanic Cloud (LMC) (jSreekumar et al.lll992l ). The LMC is an ex- 



cellent target for studying the link between cosmic-ray acceleration and gamma- 
ray emission since this galaxy is nearby (bringing the sources fluxes in reach 
of modern gamma-ray telescopes) and since the system is nearly seen face-on 
(avoiding the superposition of sources along the line of sight that hampers stud- 
ies in our own Galaxy). In addition, the LMC is rather active, housing many 
supernova remnants, bubbles and super bubbles and massive star forming re- 



gions that are all potential sites of cosmic-ra y acceleration (jBiermannI |200J; 
Cesarskv &: Montmerld[l983l : iBinns et al.ll2007l ). 



The Large Area Telescope (LAT) aboard the Fermi Gamma-Ray Space Tele- 
scope (FGST) provides now the capabilities to study diffuse gamma-ray emission 
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from nearby galaxies in dept h, and of the LMC in particular ( Digel et al. 20001 : 
Weidenspointner et al. 20071 ). We report here on the initial analysis of observa- 



tions taken in the course of the first year's all-sky survey by the LAT 



2. Observations 

The LAT is the primary instrument on the FGST satellite which has been 
launched from Cape Canaveral on June 11th, 2008. The LAT is an imaging, 
wide field-of-view, high-energy gamma-ray telescope, covering the energy range 
from below 20 MeV to more than 300 GeV. The LAT is a pair-conversion tele- 
scope with a precision tracker made of a stack of 18 x,y silicon tracking planes 
and a calorimeter made of 96 CsI(Tl) crystals. The tracker array is covered by 
a segmented anticoincidence shield allowing for the rejection of charged particle 
backgrounds. 

The LAT has a large ~ 2.5 sr field of view, and compared to earlier gamma- 
ray missions, has a large effective area (> 7000 cm^ on axis at ~ 1 GeV for the 
event selection used in this paper), improved angular resolution (~ 0.5° 68% 
containment radius at 1 GeV) and low dead time (~ 25 //s per event). The Icr 
energy resolution in the 100 MeV - 10 GeV energy r ange is better t han ^ 10%. 
A detailed description of the instrument is given bvlAtwood et al.l (|2009l ^. The 
on-orbit instrument calibration is presented bv lAbdo et al. " (j2009l ) . 



The data used in this work covers the period August 8th 2008 - April 24th 
2009 and amounts to 211.7 days of continuous sky survey observations. During 
this period a total exposure of 2.3 x 10^" cm s^ (at 1 GeV) has been obtained 
for the LMC region. 



2.1. Data preparation 

The data analysis presented in this paper has been performed using the 
ScienceTools version v9rll and the instrument response functions P6_V3. We 
collected all data obtained within the period August 8th 2008 - April 24th 2009 
and applied the dijfuse event class filter that has been designed to minimize con- 
tamination by instrumental backgroun d while retaining a su bstantial fraction of 
the signal. As has been pointed out by lAtwood et al.l ( 20091 ). any harsher event 



cut would not significantly improve the signal-to-noise ratio. 

We further excluded from the data all periods where the spacecraft has 
entered the South Atlantic Anomaly (SAA) and for which the spacecraft z-axis 
points more than 47° away from the zenith direction (the zenith direction being 
defined by the vector running from the Earth center through the spacecraft). 
While the SAA cut excludes periods of particular large instrumental background 
from the analysis, the latter cut excludes periods where the Earth enters the field 
of view. Furthermore, to minimize contamination from Earth albedo photons 
we exclude photons with zenith angles above 105° from the analysis. We further 
restricted the analysis to photon energies above 200 MeV where our current 
knowledge of the instrument response implies systematic uncertainties that are 
smaller than ~ 10% and where the redistribution of photons in energy due to 
incomplete energy measurements becomes negligible. 
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Figure 1. Preliminary adaptively smoothed (s.n.r. — 10) Fermi/hAT 
counts map of a 10° x 10° large region centered on the LMC for the energy 
ra nge 200 MeV - lOOGeV (greyscale). The contours show the extinction map 
of lSchlegel et al.l ()1998[ ) as tracer of the total gas column density in the LMC. 
Ten linearly spaced contour levels are plotted. The diamond in the north-east 
of the image design ates the location of the blazar CRATES J060106-703606 
(jHealev et al.ll2007l) that contributes at a low level to the gamma-ray emission 
in this area. 



2.2. Morphology 

To illustrate the distribution of observed gamma-ray photons in the LMC re- 
gion we show in Fig. [1] a counts-map of the area. The arrival directions of 
observed photons in the 200 MeV - 100 GeV energy range have been binned 
into 3' X 3' large pixels covering a 10° x 10° large area around the position 
{I, b) = (279.5°, —33.0°). The binned map has t hen been smoothed using a 2D 
adaptive Gaussian kernel smoothing technique ( Ebeling et al. 20061 ) to remove 
Poissonian noise that arises from the relatively small number of counts that have 
been registered. The signal-to-noise ratio (s.n.r.) has been set to 10 to reduce 
statistical noise variations to below ^ 10% in the image. 

We overlay as con tours on the Fermi/LAT counts map the extinction map 
of Schlegel et al. ( 19981 ) as tracer of the total gas column density in the LMC. To 
first order the extinction scales linearly with total gas column density, and we 
chose 10 linearly spaced contours to allow the reader to visually appreciate the 
distribution of gas column densities in the LMC. Obviously, a substantial fraction 
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Figure 2. Preliminary longitude (top) and latitude (bottom) photon in- 
tensity profiles of the LMC region for the energy range 200 MeV - 100 GeV. 
The solid line indicates the expected contributions from diffuse galactic emis- 
sion, diffuse extragalactic emission, instrumental background and the blazar 
CRATES J060106-703606 at (lb) = (281.04°, -29.63°) in this area of the sky. 



of the gas is found in a small area in the north of the LMC, at roughly {l,b) ^ 
(279.5°, —31.5°), which coincides with the 30 Doradus star forming region. 

The high-energy gamma-ray photons that are observed from the LMC also 
peak in this area. The photon intensity in the 30 Dor region exceeds ~ 300 
counts/deg^ while in most of the remaining regions of the LMC is remains below 
~ 120 counts/deg^ (the background rate around the LMC is around ~ 50 — 70 
counts/deg^). The excess near 30 Dor is also clearly seen in the longitude and 
latitude profiles of the photon intensity observed by LAT that is shown in Fig. [2j 

Within a rectangular box covering galactic longitudes 274° < / < 284° and 
galactic latitudes —36° < b < —30° we find a total number of ~ 1800 counts 
within the energy range 200 MeV - 100 GeV above the expected contributions 
from galactic diffuse emission, extragalactic diffuse emission, instrumental back- 
ground, and the blazar CRATES J060106-703606. These background contri- 
butions have been estimated by fitting spatial and spectral templates of their 
emission components together with a spatial template for the LMC emission 
to the data. Galactic diffuse emission ha s been mode lled spatially and spec- 
troscopically using the GALPROP model (|Stronell2007l ) version 54_59Xvarh7^| 
while for the combination of extragalactic diffuse emission and residual instru- 
mental background we assume an isotropic emission with a power law spec- 
tral distribution. CRATES J060106-703606 is modeled as a point source at 
{I, b) = (281.04°, -29.63°) with a power-law spectral distribution. For the LMC 



^ Available from the website: 



http: / /galprop. stanford.edu 
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we use the extinction map of ISchlegel et aL ( 19981 ) as spatial template from 
which we subtract a pedestal level of 0.07™ from all pixels and for which we 
set all pixels outside a radius of 4° around {l,b) = (279.65°, —33.34°) to zero in 
order to extract the LMC emission. As spectral model we assume a power law 
for the LMC. 

To describe the morphology of the high-energy gamma-ray emission from 
the LMC we first fit a point source with free position and fiux on top of the back- 
ground mo defl to our data. This results in a best-fitting point-source position of 
(/, b) = (279.58°, -31.72°) with a statistical 95% confidence error radius of 0.09° 
(the systematic position uncertainty is estimated to less than 0.02°). We note 
that this position is close to that of R 136, the central star cluster of 30 Dor, 
which is located at {l,b) = (279.47°, —31.67°), i.e. at an angular distance of 
0.11° from our best-fitted point-source location. 

The detection significance of the LMC can be estimated using the so-called 
Test Statistics (TS) which is defined as twice the difference between the log- 
likelihood Li that is obtained by fitting the LMC model on top of the background 
model to the data, and the log-likelihood Lq that is obtained by fitting the 
background model only, i.e. TS = 2(Li — Lq). Under the hypothesis that our 
model satisfactorily explains the Fermi/LAT data, TS follows a Xp distribution 
with p degrees of fr eedom, where p is the number of free parameters in the LMC 
model ( CashI [1979I ) . In the particular case of a point source with free position, 
fiux and spectral index we have p = 4 and the measured TS of 869.1 corresponds 
to a significance of 29. So". 

As next step we replace the point source model by an extended source 
model which we implement as axisymmetric 2D Gaussian shape with variable 
angular size a. In addition to the size we again fitted the position, fiux and 
power law spectral index of the source. This results in a best-fitting source 
position of (/,6) = (279.5°, -32.2°) (with a 95% confidence radius of 0.1°), and 
source extent of o" = 1.0 it 0.1°. The TS amounts to 1088.5 which is larger 
by 219.4 than the value obtained for the point-source model. Since we added 
one additional parameter (the source extent a) with respect to the point-source 
model we obtain the significance of the source extension from the Xi distribution 
to 14.80-. 

Alternatively to the geometrical models we also compare the Fermi/LAT 
data to various spatial templates that trace the interstellar matter distributions 
in the LMC. For ne utral hydro gen (H I) we use the aperture synthesis and 
multibeam data that iKim et all (j2005. ) have combined from ATCA and Parkes 
observations. For molecular hydrogen we use CO obseryatioii s of the LMC 
obtained with the NANT EN telescope ( Yamaguchi et al.ll200ll ). We further 



used the extinction map of ISchlegel et al.l ( 19981 ) (SFD) as tracer of the total gas 
column density and compare also our data to the 1 00 p,m IRIS map that has been 
obtained by reprocessing the IRAS survey data (jMiville-Deschenes fc Lagache 



^From now on we call the combination of the GALPROP model, the isotropic model and 
the CRATES J060106-703606 point source the background model of our analysis. The free 
parameters of this background model are the normalization of the GALPROP model, the 
intensity and spectral slope of the isotropic component, and the flux and spectral slope of the 
CRATES J060106- 703606 point source. 
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Table 1. Comparison of maximum likelihood model fitting results (see text 
for a description of the models). Column 1 gives the model used to the fit the 
LMC data, column 2 gives the TS value of the fit, and column 3 specifies the 
number of free parameters of the LMC model. 



LMC model TS Parameters 



Point source 


869.1 


4 


2D Gaussian source 


1088.5 


5 


H I gas map 


1173.4 


2 


CO gas map 


932.2 


2 


H I + CO gas maps 


1176.1 


4 


SFD extinction map 


1179.6 


2 


IRIS 100 fim infrared map 


1179.1 


2 



20051 ). The results of this comparison are summarized together with that of the 
geometrical models in Tabled) 

The best fits are obtained for the SFD extinction map and the IRAS 100 /im 
infrared map which give TS values of 1179.6 and 1179.1, respectively. For 2 free 
parameters (the total flux in the map and the spectral index) this corresponds 
to a detection significance of 34. 5o". An almost equally good fit is obtained 
using the neutral hydrogen map. Fitting instead the CO map to the LAT data 
provides a rather poor fit, suggesting that the gamma-ray morphology differs 
from that of molecular gas in the LMC. Fitting the H I and CO maps together 
to the data confirms this result since the fit attributes 97% of the total flux to 
the H I component. Correspondingly, the TS increase with respect to fitting the 
H I gas map alone is also negligible. 

The H I/SFD/IRIS 100 fim maps fit the data considerably better than 
a single point source, adding further evidence that the observed high-energy 
gamma-ray emission is extended in nature. Furthermore, the 2D Gaussian source 
model cannot reach the fit that is obtained by those tracer maps, suggesting that 
the emission morphology is more complex than a single Gaussian shape. 



2.3. Spectrum 

Using the extinction map of ISchleeel etaP (|l998l ) (i.e. our best fitting spatial 
template of the high-energy emission) we extract a spectrum of the LMC by 
fitting the data in 12 logarithmically-spaced energy bins covering the energy 
range 200 MeV - 20 GeV. Above 20 GeV, photons from the LMC become too 
sparse in our actual data set to allow for meaningful spectral points to be derived. 

Figure [3] shows the LMC spectrum that has been obtained by this method. 
Our analysis indicates a spectral steepening of the emission with increasing en- 
ergy, suggesting that a simple power law is an inadequate description of the 
data. We confirm this trend by fitting the data using a broken power law in- 
stead of a simple power law. This results in an improvement of the TS by 10.1 
corresponding to a significance of 2.7c7 {p = 2) of the spectral steepening. Fit- 
ting alternatively an exponentially cutoff power law improves the TS by 7.8 with 
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Figure 3. P reliminary spec trum of the LMC obtained by fitting the ex- 
tinction map of lSchlegel et al.l ([2998) in 12 logarithmically-spaced energy bins 
covering the energy range 200 MeV - 20 GeV to the Fermi/LAT data. Errors 
are statistical only. 



respect to the simple power law, corresponding to a significance of 2.8a {p = 1) 
of the spectral cutoff. Integrating the broken power law or the exponentially 
cutoff power law model over the energy range 100 MeV - 500 GeV gives iden- 
tical photon fluxes of (3.1 ± 0.2) x 10"^ ph cm s and an energy fluxes of 
(2.0 lb 0.1) X 10"^'^ erg cm^^s""*^ for the LMC. The systematic uncertainty in 
these flux measurements amounts to ~ 10%. 



Discussion and conclusions 



Sreekumar et al. ( 19921 ) reported the first detection of the LMC in > 100 MeV 



gamma rays based on 4 weeks of data collected with the EGRET telescope 
aboard CGRO. Due to EGRET's limited angular resolution and the weak emis- 
sion detected from the LMC, details of the spatial structure of the galaxy could 
not have been resolved. However, it has been obvious from EGRET data that 
the LMC was an extended gamma-ray source. 

Fermi/LAT allows now for the first time to clearly resolve the gamma-ray 
emission of the LMC and to attribute the emission maximum to the 30 Dor star 
forming region. While this coincidence could be taken as a hint for an enhanced 
cosmic-ray density in 30 Dor with respect to the rest of the galaxy, we note 
that a substantial fraction of the interstellar gas of the LMC is also confined to 
the 30 Dor area. Consequently, the target density for cosmic-ray interactions 
is greatly enhanced in this region which implies a corresponding enhancement 
of the gamma-ray luminosity. Whether the data do also support an enhanced 
cosmic-ray density in 30 Dor with respect to the rest of the galaxy needs a more 
detailed analysis of the observations. 
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The rather poor fit of the CO map to the LAT data suggests that the overall 
distribution of gamma-ray emission differs from that of molecular hydrogen. 
The distribution of neutral hydrogen fits the data considerably better and the 
combined fit of H I and CO maps indicates that any contribution to the gamma- 
ray emission that is correlated to the molecular gas is at best marginal. This 
agrees well with expectations since th e gas budget of th e LMC is largely (90- 
95%) dominated by neutral hydrogen ( Fukui et al.lll999l ). Consequently we are 
presently unable to determine the CO-to-H2 conversion factor, Xco? from our 
LM C data. 

Fichtel et al. ( 199ll ) performed a detailed modelling of the cosmic-ray dis- 



tribution in the LMC and predicted an integra ted > 100 MeV ph o ton f lux of 
(2.3 ± 0.4) X 10"^ ph cm-2s-i for the galaxy. iPavlidou fc Fields! (|200lh pre- 
dicted an integrated > 100 MeV photon flux of 1.1 x 10~^ ph cm~^s~^ based on 
estimates of the LMC supernova rate and total gas densities. Our observed flux 
of (3.1 lb 0.2) X 10~^ ph cm~^s~^ falls at the hide side of these estimates, yet 
given the uncertainties in the models the agreement can be judged satisfactorily. 

Further studies of the LMC with Fermi/LAT will now concentrate on the 
spectral analysis of the data, with particular emphasize on variations of the 
spectral shape throughout the galaxy. Thanks to the excellent sensitivity and 
angular resolution of the LAT, this is the first time that such studies become 
possible. And other nearby galaxies await their detection, such as the Small 
Magellanic Cloud or the Andromeda Galaxy (M31). Both should be in reach 
of Fermi and the comparative study of their diffuse gamma-ray emission should 
help to understand the impact of the environment and metallicity on the physics 
of cosmic rays. 
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